Membrane Heterogeneity: Manifestation of a Curvature-Induced Microemulsion 
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To explain the appearance of heterogeneities in the plasma membrane, I propose a hypothesis 
which begins with the observation that fluctuations in the membrane curvature are coupled to the 
difference between compositions in one leaf and the other. Because of this coupling, the most easily 
excited fluctuations can occur at non-zero wavenumbers. When the coupling is sufficiently strong, 
it is well-known that it leads to microphase separation and modulated phases. I note that when the 
coupling is less strong, the tendency towards modulation remains manifest in a liquid phase that 
exhibits transient structure of a characteristic size; that is, it is a microemulsion. The characteristic 
size of the fluctuating domains is estimated to be on the order of 100 nm, and experiments to verify 
this hypothesis are proposed. 



I. INTRODUCTION 

Certainly one of the most interesting models of the 
plasma membrane is that, rather than being homoge- 
neous, it is characterized by aggregates of saturated lipids 
and cholesterol which float, like rafts, in a sea of unsat- 
urated lipids^. An impressive array of experiments sup- 
port this hypothesis^, and limit the size of such aggre- 
gates in mammalian cells to the order of tens or hundreds 
of nanometers^"—. Experiments also limit the lifetime of 
the aggregates so that they are more readily described 
as dynamic domains^. The hypothesis remains contro- 
versial, however, due in part to a lack of a firm physical 
basis for the appearance of such domains. 

A few explanations have been put forth. One arises 
from the fact that model membrane mixtures of choles- 
terol and saturated and unsaturated lipids readily un- 
dergo separation into two liquid phases, one rich in the 
first two components, the other rich in the third^. Hence 
rafts might occur in a two-phase region and simply be a 
domain of one phase surrounded by the other. The small 
size of the domains could then be attributed to the effects 
of the cytoskeleton^ . One difficulty with this hypothesis 
is that it is known that a bilayer of composition which 
mimics the inner leaf of the plasma membrane does not 
undergo phase separation^, and that the coupling of such 
a leaf to another which does tend to phase separate pro- 
duces a bilayer in which the miscibility transition either 
occurs at a greatly reduced temperature or is eliminated 
entirelyi^iii. 

A second, related, hypothesis is that the inhomo- 
geneities occur in a one-phase region, and are simply 
those fluctuations associated with a nearby critical point 
of two-phase coexistence^^. Again the size of these fluc- 
tuations are proposed to be limited by the cytoskeltoniS^. 
This hypothesis is not only subject to the criticism that 
there may be no miscibility phase transition nearby, but 
also to the observation that fluctuations near a critical 
point exhibit little difference between their composition 
and that of the background from which they arise. Conse- 
quently they would not easily discriminate between dif- 
ferent proteins, which is the raison d'etre for the raft 
hypothesis itself. 



A third hypothesis is that the fluctuating domains are 
simply the signature of a microemulsion brought about 
by the presence of a line-active agent^i^. The difficulty 
with this proposal is that there is no obvious component 
to act as such an agent, one which would be attracted to 
the interface between the two phases. In particular it is 
clear that cholesterol is not line active as it prefers the 
phase rich in saturated lipids. Its initial addition to a sin- 
gle liquid phase brings about a phase separationiS-, that 
is, it raises the miscibility transition temperature rather 
than lowers it as a line-active agent would. As choles- 
terol is not line active, it is posited that the unsaturated 
lipids, which in biological membranes usually have one 
saturated as well as one unsaturated taili^, can in fact 
play a dual role: as a component of one of the two phases, 
and as a line-active agent between the phases^. A de- 
tailed model which encapsulates this idea and explores 
the effects on this microemulsion of the coupling be- 
tween the two leaves has been explored by Hirose et al J^. 
However such a model cannot explain the observation of 
nanoscopic domains in teranary systems that contain no 
lipids with one saturated and one unsaturated tail, as 
in the system dipalmitoylphosphatidylcholine (DPPC), 
dilauroylphosphatidylcholine (DLPC), and cholesterol^^. 

As this last example illustrates, nanoscopic domains 
can be brought about by means other than the action of 
a line-active agent. Indeed microemulsions almost invari- 
ably appear in any system which manifests modulated 
phases; they are, in general, the liquid phase to which 
modulated phases melt. It is for this reason that the re- 
cent observation of modulations of composition in giant 
unilamellar vesicles mimicking biological membranes^^ is 
so interesting; it implies that such membranes could well 
display a microemulsion. The questions that arise then, 
concern first, the nature of the interactions within the 
membrane which are responsible for the modulations, 
and second, the characteristic size of the droplets in 
the two-dimensional microemulsion to which the mod- 
ulated phases melt. A plausible scenario for the interac- 
tions which give rise to the modulated phases was pro- 
posed in the seminal work of Leibler— and of Leibler and 
Andelman^^, and I remind the reader in the next sec- 
tion of this mechanism which couples the local difference 
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in mole fractions of different lipids to the local curva- 
ture. There I also estimate the characteristic size of the 
droplets in the microemulsion expected in a bilayer with 
a cytoskeleton, and find it to be on the order of 100 nm. 
Therefore I propose that rafts can be interpreted as the 
characteristic droplets of a curvature-induced microemul- 
sion. Possible experiments to verify this hypothesis are 
proposed. 



II. THE BILAYER WITH COUPLED 
CURVATURE-COMPOSITION FLUCTUATIONS 

The proposal that fluctuations in curvature and com- 
position are coupled goes back to Leibler— and Leibler 
and Andelman3^. Their work has been extended explic- 
itly to bilayers^^"— , and I follow the last of these here. 
The basic physical idea is simple. The biological mem- 
brane consists of a plethora of distinct lipids with differ- 
ent spontaneous curvatures. Because of this variation in 
lipid architecture, fluctuations in the difference between 
local compositions in one leaf and the other couple to 
fluctuations in the curvature of the membrane. That is, 
lipids with larger headgroups and smaller tails will be at- 
tracted preferentially to the outer leaf in regions where 
the membrane bulges outward, while lipids with smaller 
heads and larger tails will be attracted to the inner leaf 
in the same regions. This affinity is directly observed in 
experiment^. 

For simplicity I consider the cholesterol and saturated 
lipid as one component in a binary system with the un- 
saturated lipid the other component. There are two order 
parameters representing the differences in mol fractions 
of these two components in the inner leaf, $i(r), and 
in the outer leaf, $o(r). It is convenient to consider the 
two linear combinations 0(r) = ($i(r) — $o(r))/2, and 
-0(r) = ($i(r) -|-$o(r))/2. The phenomenological free en- 
ergy consists of three pieces. The first is the free energy 
functional of the planar, coupled, bilayer which to second 
order in the order parameters can be written 
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The second piece is the curvature free energy, written 
here in the Monge representation in terms of /i(r), the 
height deviation from the planar configuration. 
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where A is a dimensionless coupling constant. In terms 
of the Fourier transform functions, the total free energy 
is, up to second order. 
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Within mean-field theory, one minimizes the free energy 
with respect to the membrane shape, SFfot/Sh{k) = 0, 
and substitutes the resulting height /i[0(fc)] into the free 
energy, Eq. |4l to obtain 
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This form of the free energy displays everything that is 
needed. First, the wavevector, k* , at which the compo- 
sition difference between the two leaves is softest, i.e. at 
which it shows the largest response, is the value of k that 
minimizes the coefficient of (f){k)(j)(—k). It is 
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for A > 1. 
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Thus the system is softest at a non-zero wavevector when 
A > 1. This requirement is understood as follows. The 
coupling between the curvature and the difference in com- 
positions favors a soft wavevector k* that is non-zero. 
As a structure with such a wavenumber is curved, and 
thus of larger area than when flat, this bending is op- 
posed by the surface tension 7. Further as regions of 
different compositions alternate, their occurrence is op- 
posed by a free energy per unit length between such 
regions, an energy which is proportional to the coeffi- 
cient b. Thus the curvature coupling to the composition, 
measured in terms of the competing tensions, must be 
large, i.e. A > 1. The particular consequence of this ten- 
dency to display a structure characterized by a non-zero 
wavenumber, k* ^ 0, is determined by the coefficient of 
(j){k*)4>{-k*) itself, which is equal to 
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with K the bending modulus and 7 the surface ten- 
sion. Lastly there is the coupling between the curvature, 
V^/i, and the difference in compositions between the two 
leaves, 4> 



F 



coupl 



= A(67)i/2 / d^r {y^h)4>, 



(3) 



When A is not only greater than unity, but is also greater 
than 1 + (2aK/67)i/2^ the coefficient of (j){k*)(l){-k*) is 
negative, so that the ensemble-average value of ^(fc*) is 
non-zero in equilibrium; that is the system undergoes mi- 
crophase separation. The resulting phase exhibits either 
stripes or a triangular array of domains. The possible 
occurrence of these phases was emphasized by the earlier 
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worka22i2^, and their possible manifestations in coupled 
bilayers has recently been explored^l. Indeed these struc- 
tures have been observed in some simulations of bilayers 
as predicted^^i^. Within mean-field theory, transitions 
between all phases are of first order except at a critical 
point which can occur when the average compositions of 
the two leaves are identical. However even this transi- 
tion is driven first-order by the large fluctuations in the 
directions of the wavevectors characterizing the ordered 
phases^, so that all transitions are of first order. 

The observation that I emphasize here, is that if the 
system tends toward order, but not so strongly as to 
manifest that order in microphase separation, that is if 
I + {2aK/lrfY/'^ > A > I, then the system will exist in 
a fluid phase, but a fluid that still reflects a tendency 
towards order. That tendency is manifest in its com- 
position fluctuations, which are strongest at a non-zero 
wavevector. Its structure is reflected in the structure fac- 
tor, S{k) where S{k)~^ is the coefficient of (t){k)(j){—k). 
The structure factor has a peak at k* which is non-zero 
when A — 1 > 0. The structure is also reflected in the 
correlation function, g(r), which is the inverse Fourier 
transform of S'(k). In the interesting regime in which 
(20/4/67)^/^ > A — I > 0, and for small wavenumbers 
kA;^/7 < I it is straightforward to show that g{r) be- 
haves, for large r like g{r) « r~^/^ exp(— r/^) sm{kcr+S), 
with kc and ^ given explicitly below and S a phase of 
no interest. The exponential damping, with a charac- 
teristic correlation length ^, implies that the system is 
disordered, a liquid. The oscillatory function introduces 
an additional length, fc~^, and this shows that the liquid 
is structured. It is this property of a liquid, to display 
structure at a length scale in addition to that of the cor- 
relation length, that is characteristic of a microemulsion. 
The correlation length, ^, and characteristic wavenum- 
ber, kc, are given by 
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The correlation length ^ is equal to the characteristic 
distance at the Lifshitz line at which A = I. Conse- 
quently the microemulsion structure is strongly damped. 
However the correlation length is larger than the charac- 
teristic distance for A > I which indicates that charac- 
teristic oscillations in the fluid are manifest before being 
damped out. From Eq. ^ we see that the character- 
istic distance, k~^, is on the order of, or larger, than 
(2k/7)^/^. As typical values of the bending modulus and 
the tension of a membrane in the presence of a cytoskele- 
ton are^i k « 2.7 x IQ-^^Nm and 7 w 2 x 10"^N/m, the 
characteristic size of the fluctuating regions is on the or- 
der, or greater than, 10~'^m, or 100 nm. This mean-field 
estimate indicates that the proposed mechanism could 
account for regions of the observed size. 



III. DISCUSSION 

I have proposed that inhomogeneities in the plasma 
membrane, and those observed in model membranes, are 
microemulsions brought about by the coupling of curva- 
ture to the difference in composition of the two leaves. 
This hypothesis of a microemulsion avoids the difficul- 
ties associated with ascribing such inhomogeneities ei- 
ther to phase separation or to the fluctuations associated 
with a critical point. As noted earlier, that a biologi- 
cal membrane could display a microemulsion is strongly 
indicated by the recent observation in giant unilamel- 
lar vesicles of modulations of composition in a four- 
component system consisting of distearoylphosphatidyl- 
chohne (DSPC), dioleoylphosphatidylcholine (DOPC), 
1-palmitoyl 2-oleoylphosphatidylcholine (POPC), and 
cholesterol'^. That curvature is strongly indicated 
as the mechanism responsible for bringing about the 
modulations is also evidenced by other results of this 
four-component system. With fixed mole fractions of 
DSPC and cholesterol, the relative mole fraction p = 
[DOPC] /( [DOPC] -f [POPC]) was varied. One expects 
that the difference in spontaneous curvature between 
DSPC and DOPC is greater than that between DSPC 
and POPC. Therefore increasing the fraction p from 
small values should drive the system towards a modu- 
lated phase, and this is indeed what is observed. Simi- 
larly decreasing the value of p from the modulated phase 
is expected to cause it to become unstable to a fluid 
phase, one which would appear uniform to fluorescence 
microscopy. Again this is what is observed. 

Additional support for the mechanism proposed here 
would be provided by an estimate of the value of A which 
characterizes the strength of the coupling \{b^Y/'^ be- 
tween the curvature and the difference of lipid mol frac- 
tions. It must be on the order of unity for a microemul- 
sion to occur. Leibler and Andelman^^ in their original 
paper, and later Liu et alf^S,, reasonably assume that the 
energy per unit length, A(67)^/^, should be set equal to 
k5H where 5H is the difference in spontaneous curva- 
tures of cholesterol-rich raft domains and the phospho- 
lipid background. With the coupling h of Eq. ([T]) of 
order ksT, this yields a simple expression for the di- 
mensionless coupling A = [k/ZcbT] [fcBT((57?)^/7]^/^. Liu 
et al. estimated 5H to be on the order of lO^m^^ and 
took 7 = 3.1 X 10"^N/m and k = ArnksT. These val- 
ues yield an estimate for A of about 14 which would 
imply that such membranes should always display a 
modulated phase, contrary to experiment. This nega- 
tive result would support their conclusion, which they 
reached by a slightly different argument, that the cou- 
pling between curvature and concentration fluctuations 
could not explain raft phenomena. However, if one uti- 
lizes the same difference in spontaneous curvature, the 
larger surface tension 7 = 2 x 10~^N/m of Dai and 
Sheetz^ and the same reference's smaller bending mod- 
ulus K = 2.7 X IO~^^Nm= GGksT, then one obtains the 
estimate A = 0.94. This shows that it is certainly plausi- 
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ble that the couphng has the correct order of magnitude 
to bring about the existence of a microemulsion. 

The couphng that is assumed in this paper to pro- 
duce the microemulsion, one between curvature and the 
difference in the compositions between the two leaves, 
predicts that saturated lipid-rich and -poor regions in 
the two leaves are anticorrelated. This is in line with 
the results of coarse-grained simulations of the ternary 
mixture of diarachidoylphosphatidylcholine, dilinoleoyl- 
phosphatidylcholine, and cholesterol of which the first 
two components differ markedly in curvature. The sim- 
ulations show a modulated stripe phase in which the 
stripes in the two leaves are indeed anticorrelated. This 
anticorrelation has interesting consequences for the mi- 
croemulsion. For example, an area in the outer leaf rich 
in saturated lipids and cholesterol would, due to the 
damped oscillations in composition, be bordered by a 
region rich in unsaturated lipids. These areas are anti- 
correlated with regions of the inner leaf; i.e. the above 
areas in the outer leaf would face in the inner leaf a re- 
gion rich in unsaturated lipids bordered by one rich in 
saturated lipids and cholesterol. It would be most inter- 
esting, of course, to determine whether domains, either 
in the microemulsion or modulated phases, are indeed 
anticorrelated. Perhaps this could be accomplished in 
the modulated phases by tagging the lipids on the in- 
ner and outer leaves of the vesicles with different dyes. 
Another possibility would be to observe images obtained 
from the system after being subjected to freeze fracture, 
as in ordinary microemulsions^. 



The hypothesis makes additional predictions. For ex- 
ample the ternary system POPC, DSPC, and cholesterol 
does not undergo macroscopic phase separation but does 
show nanodomains'^'*. One also knows that the ternary 
system POPC, DPPC, and cholesterol does not undergo 
macroscopic phase separations^. As the difference in 
spontaneous curvature of the two cholines in the latter is 
certainly larger than in the former, one would expect nan- 
odomains to be present, and this can be ascertained by 
Forster resonance energy transfer (FRET). It would also 
be of interest to consider the ternary systems in which 
DPPC is replaced by dimyristoylphosphatidylcholine, or 
dilauroylphosphatidylcholine as these would also increase 
the curvature difference. The dependence of the char- 
acteristic wavenumber of the domains, kc Eq. Q, on 
surface tension suggests that the domain size, obtained 
experimentally by FRET, could be varied by controlling 
the tension in experiment^. 
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